Localization and advectional spreading of convective flows under parametric disorder 
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We study thermal convection in a thin horizontal porous layer heated from below in the presence of 
a parametric disorder. Parameters are time-independent, but inhomogeneous in one of the horizontal 
directions. Under frozen parametric disorder, spatially localized flow patterns appear. We focus our 
study on their localization properties and the effect of an imposed longitudinal advection on these 
properties. Our interpretation of the results of the linear theory is underpinned by a numerical 
simulation for the nonlinear problem. For weak advection leading to upstream delocalization, the 
transition from a set of localized flow patches to an almost everywhere intense "global" flow occurs 
moderately below the instability threshold of the disorderless system. Strong advection leads to 
a washing-out of the localized flows and possible scenarios of this washing-out are discussed. The 
results presented are relevant for a broad variety of active media where pattern selection occurs. 

PACS numbers: 05.40.-a, 44.25. +f, 47.54.-r, 72.15.Rn 



I. INTRODUCTION 

The interest to the effect of localization in spatially ex- 
tended systems under parametric disorder was first initi- 
ated in the context of phenomena in quantum systems; 
this eifect, Anderson localization [ij (AL), has been orig- 
inally discussed for such processes as spin diffusion or 
electron propagation in a disordered potential. Later on, 
AL was studied in various branches of physics in the con- 
text of waves propagation in randomly inhomogeneous 
acoustic (0 and refs. therein), optical media {e.g., Q 
and refs. therein), etc. AL occurs not only for a random 
parametric inhomogeneity but also for a quasi-periodic 
one {e.g., 0, Q). The common point for the listed 
works is that they deal with conservative media (or sys- 
tems) in contrast to active/dissipative ones like in the 
problems of thermal convection. 

In the effect of parametric disorder on the exci- 
tation threshold in an active medium, ID Ginzburg- 
Landau equation has been discussed, but the localiza- 
tion effects were beyond the study scope. Noteworthy, 
the localization properties of solutions to the linearized 
Ginzburg-Landau equation are identical to the one for 
the Schrodinger equation (where AL was comprehen- 
sively discussed; for instance, see [1, |^, while for 
large-scale thermal convection it is so only under specific 
conditions. To the best of the authors' knowledge, lo- 
calization (or, in an adapted formulation, excitation of 
localized modes) in the presence of a frozen parametric 
disorder in fluid dynamics problems {e.g., thermal con- 
vection) has not been reported in the literature. 

The paper is organized as follows. In Sec.|TTl we in- 
troduce the specific physical system we deal with and its 
mathematical description, and discuss some general as- 
pects of the problem. Sec. IIIII presents the study of prop- 
erties of the formal localized solutions to the linearized 
equations, discussion of the effect of an imposed longitu- 
dinal advection (pumping) on these properties, and an- 
alytical estimations for localization exponents calculated 
numerically. Then, in Scc. lIVl we underpin our results 



obtained for formal solutions to the linearized problem 
with numerical simulation of the non-linearized equa- 
tions. In Sec. 13 results on scenarios of the advectional 
washing-out of localized convective flows complete our 
discussion of the role of the advection. Sec. lVIl finishes 
the paper with conclusions and argumentation in sup- 
port of general applicability of our results. 



II. LARGE-SCALE THERMAL CONVECTION 
IN A HORIZONTAL LAYER 

A. Basic equations 

As a specific physical system, we consider a thin porous 
layer saturated with a fluid. The layer is conflned be- 
tween two impermeable horizontal plates and heated 
from below (Fig.[T]). The coordinate frame is such that 
the (x, ?/)-plane is horizontal, z = and z = h are the 
lower and upper boundaries, respectively. The bounding 
plates are nearly thermally insulating (compared to the 
layer) what implies a flxed heat flux through the layer 
boundaries. This flux Q{x,y) is time-independent but 
inhomogeneous in space; 

Q{x, y) = Qcr(l + £^q{x, y)) , 

here Qcr is the threshold value of the heat flux for the 
case of spatially uniform heating (above the threshold, 
convective flows are excited), is a characteristic value 
of local relative deviations of the heat flux from the criti- 
cal value. For thermal convection in a porous medium the 
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FIG. 1: Sketch of the system and coordinate frame 
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effects we will discuss appear also when the heating is uni- 
form but the porous matrix is weakly inhomogeneous [Tl| 
(this inhomogeneity is inevitable in real systems). Nev- 
ertheless, we assume an inhomogeneous heating in order 
to be able to extend all our results and conclusions to 
convection without a porous matrix. In this paper, we 
consider only q — q{x) and the flows homogeneous along 
the y-direction. Also, we admit a pumping of the fluid 
along the layer. 

For nearly thermally insulating boundaries and a uni- 
form heating, the first convective instability of the layer 
is long- wavelength [13, [3 (in other words, large-scale), 
what means that the horizontal scale L of the flow is large 
against the layer height h; L = e^^h 31]. For large-scale 
convection the temperature perturbations 6 = 9{x) are 
almost uniform along z and the system is governed by 
the dimensionless equation 

6 ^ {-u9 - 9^^, ~ q{x)e^ + (e^f) ^ , (1) 

where u is the x-component of the imposed advection 
(through-flow) velocity. Noteworthy, the same equation 
governs some other fluid dynamical systems [fl uBl [3] ■ 
For Eq. ([T]) the length scale is set to L and thus the di- 
mensionless layer height = e <C 1. 

Though Eq. JT]) is valid for a large-scale inhomogeneity, 
i.e., h\qx\/\q\ ^ 1, one may set such a hierarchy of small 
parameters, namely e^(/i|(7i,|/|g|)^<Cl, that, on the one 
hand, the long-wavelength approximation remains valid, 
whereas, on the other hand, the frozen inhomogeneity 
may be represented by a 5-correlated Gaussian noise: 

qix) =qo + a^), ia^)) = 0, ia^n^')) = 2DS{x^x'), 

where go is the mean departure from the instability 
threshold of the system without disorder (shortly re- 
ferred as "mean supercriticality" ) . The noise strength 
D may be set to 1 by means of the following rescal- 
ing: {x,t,q) {D-^/^x, D-^/H, D^^^q) ([H; hence- 
forth, D = 1). As soon as the layer is practically not 
inflnite, one should be subtle with the small noise limit 
which implies a shrinking of the length of the system 
owing to the rescaling of the spatial variable. 

In specific physical systems described by Eq. ([T]), prop- 
erties of various physical fields may differ, and, simulta- 
neously, for particular physical phenomena in these sys- 
tems, properties of one or another field may play a de- 
cisive role. For instance, transport of a pollutant [2l| 
is determined by the velocity field which not necessarily 
possesses the same localization properties as the temper- 
ature field. Therefore, relations between various flclds 
and field 9 deserve attention, though all the findings of 
the paper can be provided in terms of 9 regardless to a 
specific origin of Eq. ([T]) . For a porous layer we consider, 
the fluid velocity fleld of the convective flow is 

« = ^e:.-^e^, * = /(z)7/;(a;,i), (2) 

where ipi^j t) = 9x{x, t) is the stream function amplitude 
and f{z) = (SVSS/D/i^) z{h - z) The contribution 



of imposed advection u is not presented here due to its 
smallness against v {u ^ 1 vs. v ^ h~^). In spite of its 
smallness, the advection u influences the system dynam- 
ics due to its spatial properties; u provides a fluid gross 
flux through a cross-section of the layer, whereas for the 
convective flow v this gross flux is zero [s^. When the 
gross flux of a certain flow is zero, the transport (includ- 
ing the heat transfer) by this flow is essentially less effi- 
cient compared to the case of a nonzero gross flux. The 
relationship Q between the flow and the temperature 
perturbation is speciflc for a porous layer and may differ 
from the ones for other fluid dynamical systems governed 
by Eq.©. 

B. Excitation of localized modes: general aspects 

Numerical simulation reveals only the stable time- 
independent solutions to establish in the course of the 
evolution of the dynamical system ([1]) for small enough 
u. Hence, the object of our interest are the localiza- 
tion properties of nontrivial time-independent solutions. 
When some solution is localized near the point a;o: in the 
distance from a;o, 

9 / '^'^P(~^-l^ ~ ^ol)' (2:-a;o)w>0; 
\ exp(-7+|x - xol), {x - xo)u <0. 

Here 7_ and 7+ are the down- and upstream localization 
exponents, respectively; the inverse value X± = 7±^ is 
called the "localization length" . In the region of expo- 
nential tails the solutions are small and the exponents 
may be found from the linearized form of Eq. ([1]) , which 
can be once integrated in the stationary case; 

u9 + 9"' + [qo + ^{x)]9' = const = S (3) 

(the prime denotes the x-derivative) . For u ^ the sub- 
stitution 9 9 + S/u sets S to zero. 

Let us now compare our problem to the Schrodinger 
equation, a paradigmatic model for AL, and recall some 
results from the classical theory of AL, which are relevant 
for our study. 

For u — S — 0, 9x IS governed by the stationary 
Schrodinger equation with the mean supercriticality qg 
instead of the state energy and ^{x) instead of the po- 
tential. For the Schrodinger equation with a J-correlated 
potential the AL is comprehensively studied; all states 
(for any energy) are localized {e.g., [sL loL ITol| ) . 

Our problem is peculiar not only due to the advection 
u, but also in the physical interpretation and observabil- 
ity of effects related to properties of formal solutions. 
Essentially, in the linear Schrodinger equation, different 
localized solutions to the linear problem do not mutually 
interact. In our fluid dynamical system, all these modes 
do mutually interact via nonlinearity and form a whole 
stationary flow, which may be almost everywhere intense 
for a high spatial density of localized modes. The role of 
nonlinearity for the the Schrodinger equation is addressed 
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FIG. 2: Spatial density of centers of flow excitation in sys- 
tem (Ul at M = 0. Squares: results of numerical simulation, 
solid line: dependence Q. 

in the literature {e.g., [l3ilI3) where it is reported to be 
able to lead to destruction of AL and other interesting 
effects. However, the physical meaning of the quantum 
wave function imposes strong limitations on the nonlin- 
earity form (conservation of particles, etc.), whereas in 
a fluid dynamical system similar limitations are not ap- 
plied. Therefore, the results of studies for nonlinear ef- 
fects similar to [10, II4I niay not (not always) be directly 
extended to the effects one faces in fluid dynamics. 

As additionally follows from the above, localized flows 
can be observed for a low density of excited modes. This 
is the case of large negative deviation qo of the mean 
heat flux from the critical value. To make it clear let us 
consider a local mean value 

2 i-x+l/2 

= 7 q{xi)dxi . 

' Jx-l/2 

When at the point x the value of qi (x) is positive for large 
enough / ("large enough means Z ~ 1), one may expect 
a convective flow probably to arise in the vicinity of this 
point, while in the domain of negative qi the convective 
flow is averagely damped. This criteria is rather qualita- 
tive and the speciflc value of / is not imposed by circum- 
stances; for the sake of dcflnitcness we choose I = tt [33| 
for representation of q{x) in flgures. For large negative 
qo, positive values of qi are of low probability and, there- 
fore, the domains of an intense convective flow may be 
confidently expected to be sparse. 

The last expectation is validated by numerical results 
[3^ on the spatial density v of the domains of flow ex- 
citation presented in Fig.O As qi is a Gaussian random 
number with mean go and variance 2/1, the probabil- 
ity P{qi > 0) is the error function of qo^/l/2; for large 
negative qo, P{qi > 0) w (-(7o)"X7r/)-i/2 exp(-qgZ/4). 
Fortunately, the dependence 

z^app = 0.25P(g/=i.95>0) (4) 

fairly flts the numerical results in Fig. [51 Note, the quan- 
tity u makes sense only when excitation centers are sepa- 




FIG. 3: The upper figure shows extinction of the excitation 
centers as qo decreases for a sample ^(2;); the squared stream 
function amplitude il)^{x,q(,) of the establishing pattern is 
plotted. The lower figures present fiow stream lines for qo 
indicated in the plot (the scales of x and z are different). For 
qo < —2, no flows occur as the size of the domain presented 
becomes comparable to the mean distance between excitation 
centers (cf. Fig.[2ll. 

rable, i.e., the characteristic inter-center distance is large 
compared to the localization length. In Fig. [3] we show 
typical flow patterns and demonstrate extinction of the 
excitation centers as qo decreases for a sample £,{x). As 
excitation centers becomes sparse, separate localized flow 
patterns can be apparently discriminated and considered; 
Fig.[3] suggests go ^ —1.5 to be the case. 

III. LOCALIZATION EXPONENTS 

A. Spatial Lyapunov exponents 

Let us consider the stochastic system 

e'^^, i;' = <j>, ^' = -[qo + ~ u0 , (5) 

which is equivalent to Eq. ^ with S — 0. One may treat 
the system ([5]) as a dynamic one with the spatial coordi- 
nate X instead of time and evaluate (spatial) Lyapunov 
exponents (LE) which yield eventually the localization 
exponents {e.g., in review fld\ the LE is employed as a 
localization exponent in order to estimate the localization 
length in classical AL). 

The spectrum of LEs consist of three elements: 71 > 
72 > 73- System ([5]) is statistically invariant under the 
transformation 

{u, X, 9, -0, 0) -> {-u, -X, 0, -i/j, 0) 

[recall, ^{x) possesses spatially uniform and isotropic sta- 
tistical properties and an even distribution] . Changing u 
to —u, this transformation simultaneously turns 71 (u). 
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72 (m), 73(1*) into -71 (w), -72 (u), -73("); after arrange- 
ment, one finds — 73(u) > — 72(w) > — 7i(u) which should 
be the spectrum 71 (—u) > 72 (—w) > 73 (—u)- Therefore, 

7i(9o,u) = -73(90, "u), (6) 

72(90, u) = -72 (go,-") • (7) 

Additionally, the last relation means 72(50, w = 0) = 0. 
For u = the system admits the homogeneous solution 
{6,ilj,(j)} = {1,0,0}, which corresponds to the very LE 
72 — 0. As the divergence of the phase flow of system ^ 
is zero, 71 + 72 + 73 = 0; therefore, 

72(90, u) = -7i(go,ti) +71 (go,-") ■ (8) 

The properties ([6])-([71) are demonstrated in Fig.|4^ with 
the spectrum of LEs for go — -1- Thus, due to © and 
(HI), it is enough to calculate the greatest LE 71. 

Fig. [4] shows the dependence of 71 (b) and 72 (c) on 
advection velocity u and mean supercriticality qq. For 
any u, decrease of qq leads to growth of 71, i.e., makes 
the localization more pronounced. For qq < 0, where 
localized flow patterns may occur, and non-large u, 71 
decreases as u increases. In the following we will disclose 
the importance of the fact, that 72 is positive for positive 
u. 

Let us now discuss relationships between the LEs eval- 
uated and the solution properties. When, for u — S — 0, 
some flow is localized near xq, in the distance from xq, 

f 62,- + ei(a;)e''i(^-^"), x < xq, 
9{x) w < , ^ 

i e2,+ + e3(a;)eT^("-"''), x > xq, 

where Oi(a;) and 63(3;) are bounded functions (i.e., nei- 
ther decay nor grow averagely over large distances) de- 
termined by a specific realization of noise [36j . Con- 
sidering the solution in between of two excitation do- 
mains near xi and X2 > Xi, one can combine the 
above asymptotic laws for x < xq (assuming xq — X2), 
i.e., 82,- + Qi{x) e''^^^^'^'^\ and for x > xq (assuming 
2^0 = a^i), ©2.+ + ©3(2;) e^-''^^~^i\ and obtain 

eixi<x<X2) « e3(a;)e^='(^-^i) + 63 + ei(a;)e^i(^-^^l 

(9) 

Here the amplitude of the corresponding stream func- 
tion V'(x) « 5'3(a;)e''^(=""^i)+«'i(a;)e'>'i(''-^^) [^1,3(2^) are 
bounded like ©1^3(2;)] is localized near xi and X2 with 
the exponent 71 [73(90, " = 0) = —71(90, u = 0)]. Mean- 
while, the temperature perturbation ^ is not localized 
because of constant 02 which is generally different be- 
tween different neighboring excitation regions. 

For u > (the case of u < is similar and does not 
require special discussion), the shift of the temperature 
6^6 + S/u eliminates the heat flux S. From the claim 
5 = for w 7^ 0, it follows that in the domain of the flow 
damping, where '0^0, the temperature perturbation 9 
tends to zero as well. Indeed, solution ([9]) takes the form 

e{xi<x<X2) ~ Ssix) e^^(^-=^i) 

+ 62(2;) e'^^f^-^^) -I- ei(x) e^i(^-^^) (10) 
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FIG. 4: Spectrum of Lyapunov exponents of system ([S} is 
numerically evaluated by a direct simulation. Spectrum 71 > 
72 > 73 at go = —1 demonstrates the spectral properties in 
plot (a). The dependencies of 71 and 72 on u and qo are 
presented in plots (b) and (c), respectively. 

[as 72(90 < 0, u > 0) > 0, the mode corresponding to 72 
is localized near X2], and, for u ^ 0, the contribution 
of the second term tends to in the distance from exci- 
tation domains. Hence, in the presence of the advection 
the temperature perturbations are localized as well as the 
fluid currents. On the other hand, now the mode corre- 
sponding to 72 makes a nonzero contribution to the flow: 
V'(x) « 'i'3(x)eT^(^-^i)+^'2(a;)e^^(^-^2)-H«'i(a;)e''i(^~"=^). 

As a result, for 72 > (it > 0), the flow is localized 
down the advection stream (the right flank of the local- 
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FIG. 5: (a): Greatest growth exponent /i of the mean-square 
temperature is plotted for system ([5]). (b): Exponent fi 
(squares) is compared to doubled greatest Lyapunov expo- 
nents 271 of system ([5} (solid line) and of the noiseless sys- 
tem (dashed line). Though for go < the trivial solution of 
the noiseless system is absolutely stable, the latter exponent 
formally calculated for system ([5]) with ^{x) = provides a 
rough approximation to the Lyapunov exponent of system ^ 
in the presence of disorder. 



ized flow) with the exponent 73; 7_ — I73I. On the left 
flank of the excitation domain, two modes appear: 

^Pix < X2) « ^2ix) e^^("~"^) + *i(x) e^i("-"^\ (11) 

For moderate u, '^i{x) are 5*2(2;) are comparable; hence, 
the mode 'i'i{x)e'^^'--^~^^^ rapidly "disappears" against 
'i'2{x)e''^^^~^^'^ as one moves away from X2, and the 
upstream localization properties are determined by 72; 
7+ = 72. For w = 0, the function ^2(2;) = 0, and, 
for small u, 4*2 (a;) remains small by continuity. In the 
last case the flow (|lip considerably decays in the domain 
where the 71 -mode remains dominating over the small 72- 
one, and this mode determines the upstream localization 
length; 7+ = 71. 



B. Growth exponents of mean-square values 

For an analytical estimation of the greatest LE one 
may consider behavior of mean-square values of fields 
(see i; e.g., in the paper (Tsj , such an approach was uti- 
lized for approximate analytical calculation of the Lya- 
punov exponent for a stochastic system similar to sys- 
tem ([5])). Specifically, we make use of the following par- 
ticular result of which is valid for a linear system of 
ordinary differential equations with noisy coefficients: 

y'l = Lij Vj + £.{x) r,j y, . 

For normalized Gaussian (5-correlated noise S,{x) 
[{S,{x)) = 0, {i{x + x')i{x)) = 26{x')], the mean values 
{yi) (averaged over noise realizations) obey the equation 
system 



{y.,)' = {h + T%,{y,) 



(12) 



(a simple rederivation of the particular result (jl2p can be 
found in [H]). 

For y = {9,ip,4)}, matrix = 0, and noise does not 
manifest itself. This is due to the system symmetry and 
does not characterize behavior of a particular system un- 
der given realization of noise. In order to characterize 
this behavior, one have to consider behavior of mean- 
square values (cf. |3j]). For y — {6'^,9^,d(l),ilj'^ ,ip(j),(j>'^}, 
Eqs. dni) yield 

y[ = 22/2 , 

y2 = yi + y3, 

y's^ys- [qo + i{x)]y2 - uyi , 

2/4 = 2?/5 , 

y's^ye- ko + £,{x)]y4 - uy2 , 
y'e -2[(7o + ^(x)]y5 - 2uy3 . 

Hence, 
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The growth exponent /i of mean-square values is the 
greatest real (as a mean-square value cannot oscillate) 
meaningful eigenvalue of matrix A [sj- Although the 
characteristic equation of matrix A is a 6-degree poly- 
nomial of /Lt, it is just a quadratic polynomial of u, and 
one can analytically find the surface of the greatest real 
fJ.{qo, u), consisting of two sheets, in a parametric form: 
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FIG. 6: Sample establishing time-independent solutions to Eq. ([TJ without advection for go = —3.1 [realizations of q{x) are 
represented by qi^wix)]. The lateral boundaries of the domain (x = ±50) are thermally- insulating impermeable (a,c), periodic 
(b), and isothermal impermeable (d). Dashed lines: the inclination corresponding to flow decay with the exponent 71 — 1.67, 
dashdot lines; decay with the exponent /i/2 = 1.81. 



These expressions make sense for fi > 0, go ^ 2/Sfi — 
3/i^/4+ y^2/i/3. Fig. [5] presents this surface and demon- 
strates good agreement between 271 and /i (which do not 
necessarily coincide by definition). 



IV. NONLINEAR SOLUTIONS 

In this section we 

(i) underpin the findings of Sec. lIIII with the results of nu- 
merical simulation for the nonlinear problem ([T]), 

(ii) explore the role of lateral boundary conditions for a 
finite layer, and 

(iii) discuss the consequences of the revealed fact that the 
fiow patterns may be localized up the advection stream 
with exponent 72 which can be small. 

Here our ultimate goal is to check and to show, that all 
the effects and localization peculiarities found with the 
linear system can be observed in the full nonlinear prob- 
lem [thus, we validate the argument of Sec. Ill Bl that the 
localization properties can be comprehensively studied 
by means of Eq. ([3])] . 

Although the localization properties for the stationary 
linear equation ([31) at w = are well known, the be- 



havior of the fluid dynamical system ([T]) in the absence 
of advection deserves particular attention due to above 
mentioned essential differences between this system and 
the Schrodinger equation for quantum systems. 

In agreement with expectations (Sec. Ill Bj) on observ- 
ability of localized flows for qq < 0, one can see samples 
of localized solutions for u = in Fig. [SI With thermally 
insulating lateral boundaries (Fig. [5^), the localized flow 
is excited in the vicinity of the domain, where spon- 
taneously takes positive values, and decays beyond the 
excitation domain (the background noise 9"^ ~ 10~^^ is 
due to the limitation on numeric accuracy). Thermally 
insulating lateral boundary conditions for a finite obser- 
vation domain are eventually "free" ones (no external 
impositions at the boundary) and thus correspond to a 
large (compared to the domain size) distance between ex- 
citation centers in an infinite layer. For qq — —3.1, this 
distance is indeed large, approximately 3-10'^ (see Fig. [5]). 

For periodic lateral boundary conditions (Fig.[5j3), the 
finiteness of the distance between excitation centers ef- 
fectively manifests itself (this distance is the system pe- 
riod). The system Q with vanishing noise admits the 
trivial solution 9 = q^^S x which is related to a constant 
nonzero heat flux S. In the presence of noise, this mode 
turns into the time-independent solution 9' = q^^ S -{-ipi^ 
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FIG. 7; Sample establishing time-independent solutions to 
Eq. IT} for u = 0.3, go = —3.0 [q{x) is represented by qi=n{x)] 
and periodic lateral boundaries. Dashed lines: the inclination 
corresponding to flow decay with the exponent 73 = 1.70, 
dashdot lines: decay with the exponent 72 ~ 0.134. One also 
observes multistability (i.e., coexistence of states stable to 
small perturbations) for each separate localized flow pattern; 
it can change its sign. 



(tpi) — 0, which for small S obeys the equation 

^1 +[go + e(.^)]V^i = -%'^^(a;), 

[here {£,{x)ipi) =0], i.e., tpi oc S. In the presence of noise 
this solution dominates in the domain where flows are 
damped {qi < 0). In the situation presented in Fig.[BjD, 
Qq^S — [0]/Liaycr, where [6] « 2 is the temperature alti- 
tude imposed by the excited nonlinear flow, Liayor ~ 100 
is the system period (calculation domain size) . The quan- 
tity (S/qq)'^ w 4- 10~^ has exactly the order of magnitude 
corresponding to the background flow observed in the 
distance from the excitation center (Fig.[6lD). This back- 
ground flow distorts the localization and makes the local- 
ization exponents of the "tails" of the excited localized 
flow pattern hardly measurable. For a time-independent 
flow the heat flux S is constant over the layer; therefore 
the mean background flow (5') cx {6') =^ remains 
constant over the layer even for a large number of exci- 
tation centers. 

For thermally insulating lateral boundaries the heat 
flux S related to the background flow should decay to 
zero. Therefore with thermally insulating lateral bound- 
aries the time-independent localized flow patterns occur 
without the background flow for an arbitrary number of 
excitation centers (Fig.[BJ;). In its turn, with isother- 
mal lateral boundaries (specifically, maintained under 
constant temperature 9 — 0) the intensity of the back- 



ground flow decreases as the layer extends. The rea- 
son is that, for practically separated localized flow pat- 
terns, there is multistability between solutions with dif- 
ferent signs of the temperature jumps in domains of in- 
tense flows; hence, with a large number of excitation cen- 
ters, the jumps can be mutually balanced, resulting in a 
smafl value of q^^S = L'^^^^J^ji^^]] (^-ff-, in Fig.^ the 
background flow is considerably smaller than in Fig[6)D). 
Strictly spiking, for independent random [6]j, the sum 

« \/-^iaycr, and, heuce, q^^S oc L'^^i"^ . 

Sample flows for u ^ are shown in Fig.[7l In agree- 
ment with predictions of Sec. lIII Al the temperature is 
equal on the both sides of an intense flow domain, i.e., 
not only flows are localized but also temperature pertur- 
bations. And, which is more noteworthy, even for quite 
small u when 72 is also small, the localization proper- 
ties up the advection stream are determined by 72, but 
not by 71 which is valid for vanishing u (as explained in 
Sec. lIII A| . i.e., these properties may change drastically. 
For instance, in Fig. [7] for u = 0.3 the upstream localiza- 
tion length increases by factor 12 compared to the one in 
the absence of advection. The fact, that advection causes 
upstream delocalization but not downstream, is qualita- 
tively similar to the situation with barchan sand dunes, 
windward flank of which is flat [19^, 2Ql ■ 

Noteworthy, the advectionally increased localization 
length becomes comparable to the mean distance be- 
tween the excitation centers (Fig. [2]) for moderate neg- 
ative go- In this way, weak advection may lead to the 
transition from a set of localized convective flow areas 
to an almost everywhere intense "global" flow, and, for 
instance, drastically enhance transport of a nearly indif- 
fusive pollutant which is transferred rather convectively 
than by the molecular diffusion. Quantitative analysis of 
these effects is a separate and laborious physical prob- 
lem and will be considered elsewhere (this problem has 
already been initiated in [2l| but without an imposed 
advection) . 



V. ADVECTIONAL WASHING-OUT OF 
LOCALIZED CONVECTIVE FLOWS 

A. Temporal evolution of the system 

In the previous sections we discussed only time- 
independent states, whereas in the presence of a strong 
advection the establishing regimes can be expected to be 
time-dependent. Indeed, let us consider a positive spa- 
tially uniform q{x) = const > 0. In such a case a time- 
independent space-periodic pattern establishes without 
advection. Advection u turns this mode into a traveling 
wave [in the reference frame moving with velocity u the 
advectional term in Eq. ^ vanishes] . Inhomogcneities 
of q{x) destroy the space shift symmetry of Eq. ([T]) and 
affect the traveling wave regime. In our noisy case, when 
the excitation domain is very narrow, the pattern con- 
fined between the damping regions is locked and frozen. 
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while for a broader excitation domain, a successor of this 
travehng wave may occur. 

Before we could discuss results related to the tempo- 
ral evolution of patterns for irregular q{x), an impor- 
tant fact, which may seem to be counterintuitive at the 
first glance, should be addressed. In spite of disorder in 
q{x) the temporal evolution of a certain solitary local- 
ized pattern is a deterministic process but not a stochas- 
tic one [3^. In order to support this fact, let us re- 
call a situation comprehensively studied in the litera- 
ture {e.g., [13, [13), the case of a null-dimensional sys- 
tem driven by a noisy time-dependent signal. In that 
case the temporal evolution of the system is a stochas- 
tic process. Nevertheless, for an ensemble of these sys- 
tems driven by common noise the patterns in the phase 
space at a certain moment of time (so-called "snapshot 
attractors" ) either are smooth or have a fractal structure 
which is, however, proper to deterministic chaos but not 
to stochastic one, even though reference points of these 
patterns stochastically evolve in time. In this paper we 
deal with a situation which is mathematically similar to 
the recalled one up to the exchange of time and space. In 
the same fashion as, in the case of time-dependent spa- 
tially homogeneous noise, spatial patterns are regular, in 
the case of time-independent disorder in parameters, the 
temporal evolution is smooth and deterministic. 

We do not construct a rigorous proof that the tem- 
poral evolution is smooth and deterministic in our case 
(for the recalled "twin" situation a comprehensive discus- 
sion of the regularity property of snapshot attractors is 
developed in a series of works [HI). Instead, a rough in- 
tuition for this fact can be provided with the problem of 
the linear stability of the trivial state. Without disorder 
the eigenvalue spectrum of this problem is bounded from 
above (infinitely fast growing modes are not possible). 
Disorder in q{x) leads to a spreading of the spectrum. 
However, shifts of eigenvalues related to this spreading 
are finite i.e., the spectrum remains bounded from 
above. In the presence of disorder, eigenfunctions g\{x) 
corresponding to eigenvalues A become nonsmooth but 
remain finite. Thus, solution 9{x,t) = X^a 
[C\ are amplitudes] to the linear problem evolves in time 
deterministically in spite of irregularity of g\{x). The 
results of numerical simulation for Eq. ([1]) as well show 
spatially irregular patterns to evolve regularly in time. 



B. Scenarios of washing-out 

The imposed longitudinal advection u leads not only 
to the upstream spreading of localized fiows but also sup- 
presses these flows. For an everywhere intense flow it is 
not obvious (without special study) whether the advec- 
tion enhances or suppresses convection. Differently, with 
a localized flow pattern, the advection washes away the 
heat perturbations from the excitation domain and, in 
this way, suppresses convection. In the simplest case, as 
one increases the advection strength, the intensity of a 



convective flow decreases till it disappears via a pitchfork 
bifurcation (scenario I) where the trivial state becomes 
the only attractor. In scenario I, establishing flows are 
always time-independent. For some realizations of noise 
(quite often) the convective flow disappears via a Hopf 
bifurcation near the trivial state (scenario II; here the os- 
cillatory mode is a successor of the traveling wave regime 
discussed above). Moreover, in scenario II, multistability 
and hysteresis transitions as u changes are observed (the 
Hopf bifurcation is only the last stage of this complex 
scenario) . 

Let us discuss complex scenario II in details for a sam- 
ple realization of q(x) shown in Fig. llOb and periodic 
boundary conditions. Fig. [5] presents the evolution of 
the establishing temperature altitude and stability prop- 
erties of the establishing states and of the trivial one 
as u increases. First, let us note that close to point 
P the leading eigenvalue of the problem of the linear 
stability of the trivial state is small and turns complex 
(gray curves in Fig.[Hl;); therefore the system is close 
to a co-dimension 2 bifurcation where two real eigen- 
values turn to zero, but this bifurcation differs from a 
Bogdanov-Takens one [13| due to the sign inversion sym- 
metry (9^—9) of our system. Such a situation (close- 
ness to a Bogdanov-Takens bifurcation affected by the 
sign inversion symmetry of the flclds) is wide spread in 
fluid dynamics and, e.g., Knobloch and Proctor have 
encountered it for the problem of convection in an im- 
posed vertical magnetic field [1^. As on some stages 
of the bifurcation scenario the central manifold of the 
system is 2D (cf. [11|), we project system trajectories 
onto the 2D plane (||a;6'||, ||x^6'||) (henceforth, the nota- 
tion 11. ..11 = f, ...dx stands for integration over the 

calculation domain) for presentation in Fig. [5] [these lin- 
ear in 9 integral quantifiers of the pattern 9{x) are the 
simplest independent ones]. 

In the absence of advection, only time-independent 
fiows can establish. For small u, this is still valid by 
continuity (Fig. [5^). As advection becomes stronger, a 
pair of finite-amplitude oscillatory fiows (Fig.[5]D) ap- 
pears at u — ut ~ 0.754 via a tangential bifurcation 
(point T in Fig. [5]); the more intense periodic fiow is lo- 
cally stable (Fig.[9lD, sofid loop: stable hmit cycle, dashed 
curve: saddle cycle). At ush « 0.75956 (Sh in Fig.[5lD), 
the unstable manifold of the origin intersects the stable 
one and two homoclinic loops appear (Fig.O;; Fig. [TUB 
presents the temporal evolution of the homoclinic solu- 
tion); the union of these loops is the former unstable 
cycle. Simultaneously with these repelling loops, the lo- 
cally stable periodic fFigs.[5bfTUb) and time-independent 
fiows fFigs.[9blT0fe) exist (the homoclinic loops belong to 
the boundary between the attraction basins of these sta- 
ble fiows). For stronger advection, itgh < m < UHsubi 
the homoclinic loops turn into a pair of unstable cycles 
(Fig.[9ji) [3^, which shrink as u increases and collapse 
onto the time-independent solutions making them unsta- 
ble via a subcritical Hopf bifurcations at MH^ub ~ 0.775 
(Hsub in Fig. [5]). For wn^ub < < un^^p, the only stable 
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FIG. 8: Scenario of the washing-out of the locahzed convec- 
tive flow in the system ^ is numericaUy calculated for q{x) 
plotted in Fig. llOb ,. (a): the altitude of the temperature field 
[9] = ma,x6{x) — mm9{x) of the time- independent stable pat- 
tern (solid line) and the maximal value of the one for the 
stable oscillatory regime (dashed line), (b): zoom- in of the 
region marked by the rectangle in plot (a); gray lines present 
unstable regimes, (c); the black lines present the exponen- 
tial decay rate (solid line) and the frequency (dotted line) of 
perturbations of the time-independent pattern and the fre- 
quency of the stable oscillatory flow (dashed line); the gray 
lines present the exponential growth rate (solid) and the fre- 
quency (dotted) of perturbations of the trivial state. For the 
values of u marked by vertical lines with letters in plots (b,c), 
projections of the phase portraits are shown in Fig. [9] under 
the corresponding letters. 



regime is the oscillatory flow (Figs.lSb.fllOb). Within this 
range of u, at up ~ 0.832, the unstable time-independent 
flows disappear via a pitchfork bifurcation near the triv- 
ial state (P in Fig. [5]), turning it from a saddle into an 
unstable focus (Fig.lSf). For u > up, increase of u leads 
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FIG. 9: (Color online) Projections of trajectories of the sys- 
tem m onto the plane (||x6l||, lla;^^!!) are calculated for the val- 
ues of u indicated in the plots and q{x) presented in Fig.llOb.. 



to a shrinking of the limit cycle till the oscillatory flow 
diminishes to zero at mHbup ~ 1-062 and disappears via 
a supercritical Hopf bifurcation (Hgup in Fig. [5]). For 
u > wHsu ; Sill flows decay (Fig. [5^). Remarkably, the 
vicinity of trivial state belongs to the attraction basin of 
the time-independent flows for u < itgh and to the one of 
the oscillatory flow for ugh < u < un^^^ . 

The scenario of the washing-out depends on the spe- 
cific realization of noise; therefore no universal quanti- 
tative description may be developed. Nevertheless, one 
may roughly conclude that, for narrow excitation do- 
mains [i.e., where qi{x) has a spontaneous uplift which 
is high but narrow], the system preferably follows see- 
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FIG. 10: Temporal evolution (b,c,d) of temperature field for 
sample q{x) (a), (b): stable periodic solution at w = 0.95, 
(c,d): stable periodic solution and homoclinic one at Msh ~ 
0.75956, respectively, (e): stable time-independent state com- 
pared to homoclinic and two periodic patterns (two profiles 
are distant in time by a quarter-period) at ush [for time- 
dependent states, time is the same as in (c,d)]. One sees 
the difference in the temporal evolution of these patterns to 
not come from their topological difference which is not much 
pronounced. 




FIG. 11: (Color online) Integral characteristics of localized 
patterns at qo — —3, u = 0. (a): flow intensity \\9'i\\ = O^dx 
vs temperature altitude [6]; (b) and (c): critical advection 
strength Ucr vs [6] and H^^H, respectively. Blue squares: the 
convective flow disappears via a pitchfork bifurcation, red di- 
amonds: via a Hopf one. 



nario I, while, for broader ones, scenario II becomes more 
probable. 

Rough correlations between different integral charac- 
teristics of localized patterns without advection and the 
critical advection strength Ucr (at Wcr the localized flow 
disappears; for scenario II, Ucr = MH^up) can be ob- 
served for different realizations of noise (Fig. [Til 501). 
One may note quite strong correlation between [9] and 
1 1 6*^1 1 in Fig.fTTk. whereas correlations between critical ad- 
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vection strength Mcr and integral quantifiers of the flow 
in the absence of advection are not so well pronounced 
(Fig.[IIb,c). 

Additionally, Fig.fTT] shows washing-out scenario II to 
be the most probable one for intense convective flows and 
very rare for weak ones. This fact is in agreement with 
the results on the problem of the linear stability of the 
trivial state for sample q{x). The dependence of the lead- 
ing eigenvalue of this problem on u (gray lines in Fig.[8j;) 
possesses a similar shape for any realization of q{x), but 
point P where the eigenvalue turns complex can be either 
above or below the u-axis. In the former case, scenario II 
occurs, whereas, in the latter one, the localized flow dis- 
appears before the time-dependent regimes can be ex- 
cited, the system follows scenario I. On the other hand, 
the less intense flows (without advection) are proper to 
the lower positions of the discussed dependence for which 
the case of P lying below the u-axis, i.e., scenario I, is 
more probable, and, vice versa, the more intense flows 
the higher probability of scenario II. 

VI. CONCLUSION 

We have studied the localization phenomenon in the 
problem of thermal convection in a thin horizontal layer 
subject to random spatial inhomogeneity, including the 
effect of an imposed longitudinal advection. The study 
relies on the equations relevant to a broad variety of fluid 
dynamical systems [ill Il3, EE EB] and some other ac- 
tive media where pattern formation occurs; for instance, 
the Kuramoto-Sivashinsky equation {e.g., ^26*1 and refs. 
therein for examples of physical systems) has the same 
linear part as Eq. (jT]) with u — 0. Moreover, the ba- 
sic laws in physics are conservation ones, what quite 
often remains in the final equations having the form 
i9t [quantity] -I- V • [flux of quantity] = {e.g., in [23] the 
equation for the Eckhaus instability mode in a system rel- 
evant, for instance, to binary convection at small Lewis 
number [28*] preserves such a form) . With such conserva- 
tion laws cither for systems with the sign inversion sym- 
metry of the fields, which is quite typical in physics, or for 
description of spatial modulation of an oscillatory mode, 
the Kuramoto-Sivashinsky equation should be modified 
exactly to Eq. TO. Noteworthy, for thermal convection in 
a porous layer |11| . the frozen parametric disorder ^{x) 
may be due to random inhomogeneities of the porous 



matrix (which are inevitable in reality), while the mean 
supercriticality qq may be controlled in experiments. 

Summarizing, localized nonlinear flow patterns have 
been observed below the instability threshold of the sys- 
tem without disorder, and the dependence of the spatial 
density of the localized flow patterns on the mean de- 
parture qq from this threshold has been found numeri- 
cally [Fig. [21 approximate expression (|4])]. The up- and 
downstream localization exponents have been evaluated 
numerically (Fig. 2]) and estimated analytically [Eq. P^ ]. 
In particular, the advection has been found to result in 
the localization of the temperature perturbation in ad- 
dition to the convective flow (the former is not localized 
in the absence of advection). In agreement with theo- 
retical predictions, numeric simulation for the nonlinear 
equation has exhibited a crucial effect of the advection 
on the upstream localization properties; the localization 
length can increase by one order of magnitude for small 
finite u. Via the upstream delocalization, weak advec- 
tion may lead to the transition from a set of localized 
flow patterns to an almost everywhere intense "global" 
flow, and, e.g., essentially enhance transport of a nearly 
indiffusive pollutant (without advection the problem has 
been addressed in [21|] where a crucial role of the local- 
ized flow patterns for such a transport has been revealed 
and discussed). 

Moderate and strong advection washes out localized 
flows. This washing-out phenomenon has been studied 
for many sample realizations of noise, revealing two pos- 
sible scenarios: via a pitchfork bifurcation (I) and via a 
sequence of bifurcations finishing with a Hopf one (II). 
Scenario II has been discussed in details for a sample 
noise realization; multistability and hysteresis transitions 
were observed. 



Acknowledgments 

The authors are thankful to Prof. Arkady Pikovsky, 
Dr. Michael Zaks, Prof. Dmitry Lyubimov, and Dr. 
Sergey Shklyaev for interesting and seminal discussions 
and comments in the course of the research. The au- 
thors acknowledge the Foundation "Perm Hydrodynam- 
ics," the BRHE-program (CRDF Grant no. Y5-P-09- 
01 and MESRF Grant no. 2.2.2.3.16038), and the VW- 
Stiftung for financial support. 



[1] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 
[2] J. D. Maynard, Rev. Mod. Phys. 73, 401 (2001). 
[3] V. I. Klyatskin, Dynamics of Stochastic Systems (Else- 
vier, Amsterdam, 2005). 

[Originally in Russian: V. I. Klyatskin, Statistical de- 
scription of dynamic systems with fluctuating parameters 
(Nauka, Moscow, 1975)] 
[4] D. R. Grempel, Sh. Fishman, and R. E. Prange, Phys. 



Rev. Lett. 49, 833 (1982). 
[5] E. Borgonov and D. L. Shepelyansky Physica D 109, 24 
(1997). 

[6] J. Bourgain and W.-M. Wang, Commun. Math. Phys. 

248, 429 (2004). 
[7] M. Hammele, S. Schuler, and W. Zimmermann, Physica 

D 218, 139 (2006). 
[8] J. Frohlich and T. Spencer, Phys. Rep. 103, 9 (1984). 



12 



[9] I. M. Lifshitz, S. A. Gredeskul, and L. A. Pastur Intro- 
duction to the Theory of Disordered Systems (Wiley, New 
York, 1988). 

[10] S. A. Gredeskul and Yu. S. Kivshar, Phys. Rep. 216, 1 
(1992). 

[11] D. S. Goldobin and E. V. Shklyaeva, Large-Scale Thermal 
Convection in a Horizontal Porous Layer, Phys. Rev. E 
(submitted, 2008). [preview: arXiv:0804.2825 

[12] E. M. Sparrow, R. J. Goldstein, and V. K. Jonsson, J. 
Fluid Mech. 18, 513 (1964). 

[13] Transport Phenomena in Porous Media, edited by D. B. 
Ingham and I. Pop (Pergamon, Oxford, 1998). 

[14] E. Knobloch, Physica D 41, 450 (1990). 

[15] L. Shtilman and G. Sivashinsky, Physica D 52, 477 
(1991). 

[16] S. N. Aristov and P. G. Frick, Fluid Dyn. 24, 690 (1989). 
[17] A. S. Pikovsky and D. L. Shepelyansky, Phys. Rev. Lett. 

100, 094101 (2008). 
[18] R. Zillmer and A. Pikovsky, Phys. Rev. E 72, 056108 

(2005). 

[19] H. J. Herrmann and G. Sauermann, Physica A 283, 24 
(2000). 

[20] G. Sauermann et al, Geomorphology 36, 47 (2000). 

[21] D. S. Goldobin and E. V. Shklyaeva, Diffusion of 
a passive scalar by convective flows under paramet- 
ric disorder, J. Stat. Mech. (preprint, 2008). [preview: 
[axXi^: 0805.1518 

[22] L. Yu, E. Ott, and Q. Chen, Phys. Rev. Lett. 65, 2935 
(1990); Physica D 53, 102 (1991). 

[23] D. S. Goldobin and A. Pikovsky, Phys. Rev. E 71, 
045201(R) (2005). 

[24] Yu. A. Kuznetsov, Elements of Applied Bifurcation The- 
ory. Series: Applied Mathematical Sciences Vol. 112, 3rd 
ed. (Springer, New York, 2004). 

[25] E. Knobloch and M. R. E. Proctor, J. Fluid. Mech. 108, 
291 (1981). 

[26] D. Michelson, Physica D 19, 89 (1986). 

[27] R. B. Hoyle, Phys. Rev. E 58, 7315 (1998). 

[28] W. Schopf and W. Zimmermann, Europhys. Lett. 8, 41 
(1989); Phys. Rev. E 47, 1739 (1993). 

[29] D. S. Goldobin and D. V. Lyubimov, JETP 104, 830 
(2007). 

[30] L. P. Shilnikov, Soviet Math. Dokl. 7, 1155 (1966). 

[31] Remarkably, for a large departure from the instability 
threshold the large-scale approximation can be utilized as 
well, and the equations (in [29| they have been derived for 
an even more general problem) are simpler than the ones 
relevant near the threshold. But those equations admit 
discontinuities of spatial derivatives of the temperature 
field, and these discontinuities cannot be studied within 
the framework of the large-scale approximation. 

[32] Performing the rescaling, one should keep in mind, that 
the statistical effect of a delta-correlated noise is deter- 
mined by the integral of its correlation function which 
is influenced by coordinate stretching. In details, under 
the rescaling q = Qq (and thus ^ = QS,), x = Xx, 
and t = Ti, one flnds {^(x') S,(x' + x)) = 2DS(x) = 
2DX-^S(x) = Q'^{i(x)i(x' -f i)). In order to obtain 



a normalized noise, i.e., {^(x') ^(x' -\- x)) = 25(x), one 
should take — D/X; to preserve Eq. ((T} unchanged, 
one has to claim T = X'^ and Q = X~^. These conditions 
yield X = L>-l/^ T = D-^^^, and Q = D^/s. 
[33] In order to proof this property of v, let us consider the 
path r on the (x, 2:)-plane, which connects the point Fi 
to the point r2 on the lower and upper boundaries of 
the layer, respectively. The gross flux through this path 
is J = Jj^v ■ ndl, where n is orthogonal to P and di- 
rected into the same region as the gross flux, / is the 

length along P. Then J = J (^^^^+^'^ (~^) ) ^ 

yAaidi+a^dij^' = dT^' = - 

$(ri) — because = on the layer boundary [cf. 
Eq.@]. 

[34] The reason for such a choice is the fact, that in the noise- 
less case with all the coefficients equal 1 the fundamental 
mode is sin a;; for this mode the smallest space domain is 
of length TT. 

[35] All the numerical simulations for Eq. ^ are performed 
by means of the finite difference method (the first x- 
derivatives are central; the noiseless a;-derivatives have 
accuracy dx^; the time step dt = da:;*/ll; the noise ^(x) 
is generated in the middle between the mesh nodes). The 
space step dx = 0.1 appears to be enough small for cal- 
culations what is ensured (i) by the results in Figs. 6,7 
where the localization properties of calculated patterns 
are in agreement with the independent and much more 
accurate simulation for null-dimensional system ((5)l and 
(ii) by the fact, that statistical properties do not change 
when step da; is halved. This is due to the fact, that 
though a 5-correlated noise reduces the strong order of 
accuracy of numerical schemes, the weak order, that de- 
termines accuracy for statistical properties, is generally 
higher than the strong one. 

[36] In order to write down these expressions on the both 
sides of point xq, we take (in a general form) superposi- 
tions of all the eigenmodes corresponding to LEs 71 > 0, 
72 = (with spatially homogeneous eigenmodes 62, ±), 
and 73 < and exclude the ones that infinitely growth 
as we move away from xq. 

[37] A formal eigenvector is meaningful, when it meets the 
conditions (6^) > 0, {^^) > 0, (9!)^) > 0, > 

(e^)^ (e^)(02) > (e0)^ and (v')(0') > m'. 

[38] In the case of inappropriate initial conditions, the first 
infinitesimally small period of evolution should be ex- 
cluded (like for the conventional heat conduction equa- 
tion where inconsistencies in boundary and initial condi- 
tions are known to decay infinitely fast). 

[39] This bifurcation of homoclinic trajectories differs a bit 
from general ones described by the Shil'nikov's theo- 
rem [24. [30II due to the symmetry 9—^—6. 

[40] A precise and comprehensive investigation of the statis- 
tics on the prevalence of these scenarios and their pecu- 
liarities for different values of parameters is beyond the 
scope of the present paper 



